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Objectives 

This  research  focuses  on  multispectral  analyzing  of  weak  THz  signals  by  carbon  nanotube 
junction  arrays.  The  junction  arrays  simultaneously  determine  the  electromagnetic  field 
frequency,  polarization,  and  intensity.  Our  key  idea  involvoles  combining  of  metallic  and 
semiconducting  carbon  tubes  in  the  same  comb  array.  We  exploit  intrinsic  phase-coherent 

properties  of  carbon  nanotubes.  This  allows  improving  the  analyzer  sensitivity  up  to  ~1  fW  and 
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reducing  the  intrinsic  and  parasitic  noises  by  factor  10-10  .  The  suggested  approach  extends  the 
analyzed  spectral  diapason  from  0.5-7. 5  THz  covered  by  the  semiconducting  comb  to  7.5-50 
THz  covered  by  the  metallic  tube  comb.  Practical  outcome  of  this  work  can  be  used  for  stand  off 
detecting  of  concealed  weapon,  explosives,  toxic  substances,  and  fast  moving  projectiles. 

Practical  implementations  of  multispectral  analyzing  by  carbon  nanotube  junctions  and 
current  experimental  studies 

Sensors  of  electromagnetic  field  are  used  in  many  areas  of  science  and  technology.  These 
devices  analyze  the  spectrum  of  electromagnetic  waves  reflected  from  various  substances.  Many 
chemical  analytes,  e.g.,  explosive  or  toxic  materials  have  pronounced  spectral  features  in  the 
frequency  range  A^EF  =  0.5  —  50  THz  (see,  e.g.,  Refs.  [A-C]  and  the  references  there  in). 
Therefore  the  THz  spectral  analysis  can  be  effectively  used  for  a  stand  off  identification  of  those 
substances.  Since  the  THz  waves  (T-rays)  penetrate  through  clothing  and  paper,  the  security 
nanosensors  are  capable  of  identifying  concealed  metallic  and  ceramic  objects.  Known  medical 
implementations  of  the  T-ray  detectors  include  monitoring  the  healing  of  human  skin  and  of 
body  wounds  [4].  Most  of  T-ray  detectors  available  now  work  in  a  narrow  frequency  range,  and 
are  bulky  and  costly. 


Fig.  1.  A  comb  array  of  gated  carbon  nano  tube  junctions.  The  setup  consists  of  two 
parts:  the  comb  of  semiconducting  nanotubes  is  used  for  spectral  analyzing  of  an 
external  electromagnetic  field  in  the  diapason  0.5-7. 5  Thz  while  the  metallic 
nanotube  comb  works  in  the  upper  diapason  7.5-50  THz. 

This  motivates  an  increased  interest  to  nanoscale  devices  (nanosensors),  which  have  better 
sensitivity  and  wider  frequency  range.  The  nanosensors  are  expected  to  be  more  effective, 
robust,  cheaper  in  mass  production,  and  durable  [2,4,5,8,A-D].  The  sensing  principles  exploit  a 
general  fact  that  electron  transport  in  carbon  nanotubes  is  strongly  affected  by  external  fields  and 
environment  [1-8,A-D].  The  central  element  of  the  THz  multispectral  analyzer  is  a  comb  array  of 
semiconducting  and  metallic  carbon  nano  tubes.  The  comb  is  placed  on  a  dielectric  substrate  (in 
current  experiments  it  is  SiCE).  Each  of  the  nanotube  in  the  comb  is  controlled  individually  by 
the  bias  source  drain  and  gate  voltages.  The  bias  voltage  control  provides  an  extra  flexibility  to 
the  spectral  analyzer  setup  since  its  sensitivity  in  the  whole  spectral  range  is  optimized 
dynamically.  The  comb  array  consisted  of  two  different  types  of  carbon  nanotube  junctions: 
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semiconducting  tube  junctions  and  the  metallic  tube  junctions  in  the  Tomonaga-Luttinger  liquid 
(TLL)  state.  In  this  way  one  extends  the  analyzed  spectral  diapason  from  0.5-7. 5  THz  covered  by 
the  semiconducting  comb  to  7.5-50  THz  covered  by  the  metallic  carbon  nanotube  comb.  Besides 
our  method  exploits  intrinsic  phase-coherent  properties  of  carbon  nanotubes  which  allows 
improving  of  the  frequency  resolution,  and  reducing  of  the  intrinsic  noises  considerably. 

The  relevant  experimental  work  had  been  conducted  at  Georgetown  University  by  the 
Prof.  P.  Barbara  group.  Active  experimental  research  of  carbon  nanotube  THz  sensors  had  been 
conducted  also  in  Japan  by  Y.  Kawano  (see  Refs.  [35-37]).  In  Fig.  2  we  show  their  recent 
experimental  data  which  suggest  that  the  carbon  nanotube  junctions  work  as  sensors  of  THz 
field. 


Y.  Kawano  (2007-8) 


Fig.  2.  (a)  Experimental  data  by  Y.  Kawano  et  al.  (see  Refs.  [35-37]).  That  group  used 
semiconducting  carbon  nano  tube  junctions  serving  as  a  quantum  well  with  the  single 
electron  tunneling.  Their  devices  detect  weak  THz  fields  having  minimal  power  as  low 
as  ~1  fW  in  the  frequency  range  1-5  THz.  The  plots  above  illustrate  evolution  of  the 
single  electron  tunneling  resonances  under  the  THz  irradiation. 

The  authors  of  Refs.  [35-37]  fabricated  and  tested  the  carbon  nanotube  quantum  wells  which 
exploit  principles  the  photon-assisted  single  electron  tunneling.  Their  devices  detect  the  external 
electromagnetic  field  signals  as  weak  as  ~1  fW  in  the  frequency  diapason  1-5  THz.  Nevertheless 
the  available  spectral  range  obtained  in  works  [35-37]  is  too  narrow  while  their  working 
temperature  T<10K  is  too  low  yet  for  many  practical  purposes.  Therefore  our  work  aimed  at 
extending  of  the  carbon  nanotube  sensor  capabilities.  We  studied  possibilities  to  widen  the 
multispectral  analyzer  frequency  diapason,  to  reduce  the  noise  level,  and  to  improve  the  detector 
sensitivity. 

Our  research  findings 

1.  Sensing  an  electromagnetic  field  with  photon-assisted  Fano  resonance  in  a  two-branch 
carbon  nanotube  junction. 

A  Fano  electromagnetic  field  sensor 

In  a  typical  setup  one  measures  [2-6,8,A-D]  changes  of  electrical  characteristics  of  a  single  wall 
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carbon  nanotube  (SWCNT)  section  C  attached  to  nonnal  metal  electrodes  N  ,  and  N2  (see  the 
sketch  in  Fig.  3a  where  FSD  denotes  the  source-drain  bias  voltage  and  VG  is  the  gate  voltage), 
which  constitutes  the  nanotube  junction  (NJ).  The  NJs  of  various  configurations  are  subject  to 
extensive  experimental  and  theoretical  study  [2-6,8,A-D].  In  Fig.  3  the  nanotube  C  resides  on  the 
substrate  S  and  is  coupled  to  the  N  electrodes  via  the  Schottky  barriers  I  arising  due  to 
differences  of  electron  concentrations  in  C  and  N 1 9  .  Because  electron  motion  in  the  C  section  is 

quantized  [9],  the  carbon  nano  tube  junctions  is  quoted  as  a  one-dimensional  (ID)  quantum  dot. 
The  quantized  states  fonned  inside  the  dot  cause  various  resonances  experimentally  observed  in 
transport  characteristics  of  the  carbon  nanotube  junctions  [1,2,5,6,8,9,10]  and  quantum  dots 
[10,41].  A  steady  state  quantum  interference  between  localized  electron  states  in  C  and 
continuum  states  in  N  electrodes  results  in  sharp  Fano  peaks  and  dips  pronounced  in  the  electron 
transport  characteristics  [5, 6, 8, 9]  of  the  NJ  nanodevices. 


Fig.  3.  (a)  A  carbon  nanotube  junction  Fl§-  4-  The  stationary  transmission  coefficient  T6: 
(NJ).  (b)  A  A-sensor  of  the  EF  with  (a)  for  different  magnitude  of  the  lmewidth  T;  (b) 

electric  field  components  Ex  and  Ey.  for  vanous  field  Polarizations;  (c)  for  various  level 

energy  differences  A;  and  (d)  for  various  total  level 
energies  s. 

In  this  part  of  work  we  studied  the  influence  of  an  external  electromagnetic  field  (EF)  on 
the  Fano  resonance  in  an  carbon  nano  tube  junction.  We  will  see  that  an  external  electromagnetic 
field  induces  photon-assisted  tunneling  processes  across  the  nanotube  junction.  This  strongly 
affects  sharp  quantum  mechanical  resonances  taking  place  in  the  nanotube  junction.  The  ID 
nature  of  the  electron  transport  in  SWCNT  is  exploited  to  detennine  polarization  of  the 
electromagnetic  filed.  The  above  phenomenon  is  suggested  for  sensing  an  external  THz  signals. 
Theoretical  model 

We  consider  a  two-branch  carbon  nanotube  junction  exposed  to  an  external  EF.  The  junction  is 
composed  of  carbon  nanotube  sections  C  1(r)  coupled  to  metallic  electrodes  N 1 2  3  via  electron 

tunneling.  Quantum  interference  between  localized  electron  states  in  the  carbon  nanotube 
branches  and  continuum  states  in  the  metallic  electrodes  causes  sharp  Fano  resonances.  The 
sharp  peaks  are  affected  by  an  external  a.c.  field  which  induces  tunneling  processes  accompanied 
by  emission  and  absorption  of  photons.  We  will  see  that  the  photon-assisted  tunneling  results  in 
appearance  of  satellite  peaks  in  addition  to  the  main  Fano  resonance.  In  this  paper  we  implement 
an  illustrative  and  tractable  model,  which  describes  the  photon-assisted  phase-coherent  electron 
transport  across  a  two-branch  SWCNT  junction.  Despite  its  simplicity,  our  model  adequately 
captures  basic  physical  features  of  photon-assisted  electron  transport  under  the  conditions  of  the 
Fano  resonance.  The  SWCNT  branches  C  1(r ^  are  presumed  to  be  "clean",  which  implies  that  the 
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electron  movement  is  ballistic  (\  >  L}^  ,  where  l]  is  the  electron  elastic  mean  free  path  and 
Llir  i  is  the  length  of  the  nanotube  section).  The  C  1(r  branches  are  misaligned  by  angle  i)  and 

are  coupled  to  each  other  indirectly  via  the  common  normal  metal  electrode  N  x ,  as  shown  in  Fig. 
3b. 


Fig.  5.  A  section  of  carbon  nanotube,  its  crystal 
lattice  structure,  a  Brilloin  zone,  and  the  3d  plot  of 
electron  energy  versus  momentum. 


The  A  -shape  of  the  nanosensor  pursues  a  double  purpose:  (i)  one  measures  the  d.c.  conductivity 
G(Vsd,Vg)  in  each  of  the  C  sections  independently;  (ii)  the  ID  nature  of  electron  transport  in 

the  C  sections  allows  simultaneous  independent  measurement  of  components  Ex  and  Ey  (see 
Fig.  3b)  of  the  external  electromagnetic  field.  The  finite  misalignment  angle  $  assures  that  the 
a.c.  bias  V^'t{-h)  [here  indexes  l(r )  refer  to  the  left  (right)  branches  of  the  A -sensor,  while 
the  indexes  t(b)  are  related  to  the  top  (bottom)  Schottky  barriers]  is  induced  on  the  NJ  s 
Schottky  barriers  by  a  certain  projection  of  the  electric  field  vector  E  =  E  x  +  Eyy . 

In  the  geometry  of  Fig.  3b  one  obtains  ^  oc  Ex  sin  0  +  Ey  cos  0  and 
oc  —Ex  sin(6>  +  l)  —  Ey  cos (9  +  t) .  An  arbitrary  EF  polarization  ^  Vf^hl  gives 

different  probabilities  of  elementary  photon-assisted  tunneling  in  each  of  the  two  branches  C  ( 
and  C  r . 

The  inter-level  spacing  is  estimated  as  8e  —  hvF  /  2 L  [where  vF  is  the  Fermi  velocity 
in  the  carbon  nanotube,  vF  —  8  x  105  m/s,  and  L  is  the  nanotube  length, 
L  =  2  x  (102  —  103)  nm].  The  level  energies  £.  are  controlled  by  changing  the  charge  carrier 

concentration  when  applying  the  gate  voltages  to  C  ■ .  Since  only  one  level  e-  supports  the 
transport  across  the  i  -th  dot  (if  either  the  inter-level  spacing  8e  or  the  level  energy  difference 
A  is  larger  than  the  bias  voltages  ESD  and  VG ,  linewidth  F ,  and  temperature  T  ),  no  intradot 
Coulomb  interaction  is  present.  In  this  way  one  neglects  the  charge  accumulation  on  the 
nanotube  sections,  and  accounts  just  for  one  type  of  the  carriers.  This  is  justified  when 

max{r,e2  /  2(7}  <<  min{<5e,  A},  i.e.,  when  electron  propagation  inside  the  C-sections  is  well 
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quantized  and  the  NJ  capacitance  C  is  not  small.  Though  the  above  assumption  is  very  simple,  it 
allows  deriving  of  tractable  analytical  expressions,  which  appears  to  be  quite  convenient  for  a 
qualitative  and  quantitative  analysis. 

Photon-assisted  Fano  resonance 

For  the  sake  of  simplicity,  we  ignore  the  dependence  of  the  tunneling  matrix  elements  tsi  versus 
the  energy  variable  e  and  electron  momentum  p .  Besides,  we  assume  that  the  tunneling  matrix 
element  tsi  is  the  same  for  all  the  junctions,  i.e.,  tm  =  |f| .  Due  to  electron  scattering  on 
atomic  defects  and  tunneling  with  rates  If  and  Ys  respectively,  the  dot  energy  levels  acquire 
a  finite  linewidth  I\  +  T  s ,  where  f\  =  2tt  1 1 1  Ng  ( 0 ) ,  and  Ns  ( 0 )  i  s  the  density  of  states  in  the 
s  -th  N  lead.  Our  model  assumes  no  direct  interaction  (neither  Coulomb  repulsion,  nor  tunnel 
coupling)  between  the  sections  C  l  and  C  r  .  The  steady  quantized  states  are  coupled  to  each 
other  rather  indirectly  via  continuum  electron  states  in  the  N  leads.  We  will  see  that  a  finite 
external  electromagnetic  field  produces  series  of  sharp  satellite  singularities,  which  arise  in  the 
C'(VSD )  -curves  in  addition  to  the  steady  state  resonances.  The  position  and  magnitude  of  the 

satellite  singularities,  which  constitute  the  photon-assisted  Fano  resonance,  are  uniquely  related 
to  the  parameters  of  the  external  electromagnetic  field. 

Our  calculation  results  for  the  steady  state  transmission  coefficient  T0  ( e )  are  presented 

in  Fig.  4a-d.  In  Fig.  4a  we  plot  TQ  (e)  for  (3t  —  /3r  —  1 ,  £  —  1,  A  —  0.2,  Y  —  0.05  (curve 
1),  r  =  0.25 (curve  2),  and  Y  —  0.6 (curve  3).  One  can  see  that  Tq[e')  vanishes  at  e  =  e. 
One  may  also  notice  a  well-pronounced  double-peak  singularity  affected  by  asymmetry  of  the 
setup.  The  transmission  becomes  ideal  [i.e.,  TQ  (e)  =  1]  at  e  —  e  =  A.  Fig.  4b  illustrates  the 

effect  of  transmission  asymmetry  when  /3j  ^  (3 hr .  All  the  parameters  used  for  calculations  in 
Fig.  4b  are  the  same  as  in  the  former  Fig.  4a  except  Y  —  0.3  and  (3r  —  0.01  (curve  A), 
/3.  =  0.5  (curve  B),  and  (3r  —  1  (curve  C).  In  Fig.  4c  we  show  T0  (e)  for  various  level  energy 
differences  A  =  0.02  (curve  I),  A  =  0.2  (curve  II),  A  =  0.8  (curve  III).  All  other  parameters 
are  the  same  as  before.  In  Fig.  4d  we  plot  T{)  (e)  for  different  values  of  the  sum  level  energy 
e  =  0.5  (curve  a),  e  ==;  1  (curve  /)),  and  e  —  1.5  (curve  7).  The  partial  transmission 
coefficient  Tm  ( £ )  for  an  arbitrary  m  -photon-assisted  process  is  obtained  from  TQ  ( e )  with 
shifting  the  energy  argument  by  mMl . 

The  calculation  results  for  the  A  -sensor  conductivity  G  versus  the  source-drain  bias 
voltage  Fsd  (in  units  of  e)  are  shown  in  Fig.  6.  For  the  sake  of  simplicity,  we  assume  that  the 
a.c.  and  d.c.  biases  are  the  same  for  all  the  Schottky  barriers,  which  is  the  case  for  a  symmetric 
NJ  and  Ex  =  Ey  .  Then,  in  all  the  equations  one  sets  y^YYb)  _  anci  V^)'t<'b}  =  VSD . 

Fig.  6a  shows  the  A  -sensor  conductivity  for  eVac  =  0.1  (curve  1),  eVac  —  0.2  (curve  2), 
and  eVac  —  0.5  (curve  3).  Other  relevant  parameters  are  j3l  —  1,  f3r  =  1,  e  =  1,  A  =  0.07 

,  r  =  0.03  ,  and  Yl  —  0.4  .  One  can  see  that  the  a.c.  EF  induces  series  of  satellite  double-peak 
singularities,  with  positions  and  relative  amplitudes  depending  on  the  EF  intensity.  When  the 
intensity  of  the  external  EF  is  low,  the  photon-assisted  tunneling  signatures  are  weak:  only  a 
weak  double -peak  singularity  appears  at  the  bias  voltages  VSD  =  (s  ±  (f2  +  A))  /  e  [see  curve 
1].  Splitting  between  the  main  peak  and  the  first  adjacent  satellite  peak  in  this  case  is  Ml  /  e. 
However,  as  the  a.c.  field  amplitude  Eac  increases,  the  amplitude  of  the  main  peak  decreases 
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while  the  satellite  double-peaks  become  much  sharper  (see  curves  2  and  3).  We  emphasize  that 
here  we  neglected  the  non-equilibrium  effects  despite  the  a.c.  field  is  strong.  This  is  justified 
since  Or  <<  h  .  The  A  -sensor  conductivity  at  different  Q  is  presented  in  Fig.  5b,  where 
curve  A  corresponds  to  O  =  0.1;  curve  B  to  Cl  —  0.4,  and  curve  C  to  Cl  =  0.7.  Other 
relevant  parameters  are  (3t  =  (3r  —  1 ,  eVac  =  0.3,  e  =  1,  A  =  0.07,  and  T  =  0.03. 

Frequency  range  and  sensitivity. 

Using  the  above  data  one  may  summarize  the  EF  sensing  method  as  follows:  (z)  the 
electromagnetic  field  frequency  is  readily  extracted  from  splitting  of  the  double  peak  structure 
(Fig.  5a);  (zz)  the  field  vector  components  Er  „  (see  Fig.  3b)  are  determined  by  the  ratios  of  the 
main  and  the  satellite  double  peak  heights;  (zzz)  field  polarization  is  obtained  by  comparing  the 
satellite  peak  heights  in  the  conductivities  G,  and  G  of  the  two  adjacent  NJ  branches.  One 
utilizes  the  fact  that  an  arbitrary  EF  polarization  gives  different  sequences  of  peak  heights  in 
Gj  of  the  two  NJ  branches  l  and  r  .  From  a  measurement,  one  detennines  ratios  of  the  main 

and  the  satellite  double  peak  heights  nfi  for  each  of  the  branches  ji  —  x,y .  Then,  one  finds 

Exy  from  two  transcendental  equations  /  Jq(cx  )  =  k  ,  where  n  —  x,y .  The  x  - 

and  y  -components  of  the  a.c.  electromagnetic  field  vector  E  are  readily  obtained  from  the 
universal  dependence  a(Kxy).  The  function  a(K  )  [where  p  =  arctan(«',2.  /  Ky)]  is  given  in 
Fig.  6c. 

The  sensitivity  diagram  in  Fig.  6c  is  obtained  for  cuD  =  10,  A  =  0.3,  pF  —  1,  =  1. 

A  =  1.  One  can  see  that  the  temperature  broadening  due  to  coupling  of  electrons  to  phonons 
significantly  restricts  sensitivity  of  the  Fano  nanosensor  at  high  frequencies  (which  for  the  used 
carbon  nanotube  junction  parameters  are  below  1.7  -  17  THz)  and  temperatures  T  >  7 A 
(which  corresponds  to  80-800  K  for  the  same  device  parameters).  The  above  calculations  were 
made  under  the  assumption  that  the  external  field  amplitude  is  weak.  We  have  actually 
disregarded  non-equilibrium  effects  along  with  the  changes  in  the  electron  bandstructure  due  to 
the  photon-assisted  tunneling.  However,  the  above  assumptions  become  invalid  if  the 
electromagnetic  field  is  strong.  In  the  latter  case,  the  non-equilibrium  phenomena  influence  the 
nanosensor  characteristics  substantially.  Strong  electromagnetic  field  changes  the  electron 
distribution  function  /  ,  affecting  electron  transport  across  the  NJ  directly.  Besides,  if  the  C 

section  is  longer  than  the  electron  diffusion  length,  it  causes  instability  in  respect  to  the 
appearance  of  inhomogeneous  regions  along  C.  In  an  experiment  one  may  use  SWCNT  as  the 
left  (right)  shoulders  C1(r^  and  Pd  as  the  N  leads.  The  stationary  Fano  resonance  is  tuned  by 

applying  the  gate  voltages  VG ^  to  get  an  optimal  value  of  A .  The  gate  voltage  influence  is 
understood  when  one  considers  that  the  inter-level  energy  spacing  6e  —  hvF  /  2 L  depends  on 

the  Fenni  velocity  vF  =  (h  /  m*)(37r2n)1/3 .  The  electron  effective  mass  m*  and  the  electron 
concentration  n  both  can  be  tuned  by  applying  appropriate  gate  voltages.  The  electron  carrier 
concentration  is  altered  by  the  gate  voltage  VG  applied  to  the  nano  tube  junction  in  the  way 
[5, 6, 8, 9]  as  shown  in  Fig.  3a.  The  effective  electron  mass  is  controlled  using  the  split-gate 
configuration.  The  split  gate  voltage  USG  affects  the  electron  bandstructure  of  a  "helical" 
nanotube  already  at  E  >  10  V//z  m(l  mV/A).  Other  types  of  semiconducting  nanotubes  (e.g., 
zigzag)  require  stronger  fields  (>  0.2  V/A)  to  induce  a  metal-semiconducting  transition  when  a 
finite  energy  gap  is  opened  in  the  electron  excitation  spectrum.  For  typical  lengths 
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£1(rj  =  0.2  —  2  /_/  m  of  C1(r)  ,  one  finds  that  A  =  1  —  10  meV  (0.24-2.4  THz).  At  these 

parameters,  the  A  -sensor  works  in  the  THz  frequency  range  at  temperatures  T  =  1  —  100  K. 
A  correspondence  between  the  field  intensity  and  yTAA6)  js  obtained  from  a  self-consistent 

numeric  solution  of  Maxwell  equations,  completed  by  equations  for  electron  spectrum, 
conductance,  junction  transmission,  and  the  plot  in  Fig.  5c  addressing  the  particular  setup. 
Verification  is  made  using  calibrating  measurements.  In  this  way,  the  photon-assisted  Fano 
resonance  in  the  A  -shaped  junctions  allows  to  detennine  the  external  electromagnetic  field 
parameters. 


a)  G  |l 

hQ/e 


Fig.  6.  Calculation  results  for 
the  linear  conductivity 
G(Vsd)  of  the  A-sensor.  (a) 
The  G(Vsd)  -curves  for  three 
different  external  EF 
amplitudes  Vac-  The  splitting 
between  the  satellite  peaks  is 
oc  MU  e .  (b)  G(Vsd)  for  three 
different  field  frequencies,  (c) 
The  universal  dependence  a  = 
eVac/Q  versus  K=JC(a)/Jo  (a) 
where  Jn  is  the  Bessel 
function  of  order  n.  (d)  The 
A-nanosensor  with  antennas 
attached  to  the  NJ-branches. 


Indeed,  the  sketch  in  Fig.  3b  is  idealized.  If  the  a.c  field  frequency  is  comparable  to  c  /  L,  one 
needs  no  external  antenna,  since  one  may  utilize  the  intrinsic  antenna  effect  observed  and 
described  in  Refs.  [A]  (see  also  references  there  in).  However,  if  the  a.c.  frequencies  lay  well 
below  c  /  L,  i.e.,  when  the  EF  wavelength  becomes  much  shorter  than  the  nanotube  length,  one 
needs  an  additional  external  antenna  as  schematically  shown  in  Fig.  6d.  By  attaching  external 
antennas,  one  improves  the  reception  of  electromagnetic  signals  if  the  electromagnetic  field 
wavelength  and  the  nanotube  length  do  not  match  each  other.  The  external  antennas  may  be 
fabricated  using  lithography  methods.  Then  they  are  electrically  connected  to  the  metallic 
electrodes  of  the  A -junctions.  In  the  geometry  of  Fig.  6d,  the  a.c.  field  components  Er  and  E 

are  tuned  independently  by  each  of  the  antennas,  and  are  transmitted  to  the  corresponding  branch 
of  the  A-j unction. 

The  frequency  range  of  the  A-sensor  is  detennined  by  several  factors.  Along  with  the 
circuit  noises  and  influence  of  the  external  environment,  there  are  intrinsic  sources  of  noise. 
They  originate  from  temperature  broadening  T  of  the  Fano  resonance  peaks  due  to  the  electron- 
phonon  coupling  discussed  in  previous  subsection.  The  magnitude  of  F  depends  on  the 
temperature  and  on  the  coupling  strength  between  the  SWCNT  sections  and  metallic  electrodes. 
Another  fraction  of  broadening  is  due  to  coupling  of  the  nanotube  to  metallic  electrodes.  That 
mechanism  depends  mostly  on  the  transparency  of  the  Schottky  barrier.  The  net  broadening  is 
also  detennined  by  the  quality  of  interfaces  and  by  the  purity  of  the  electrodes. 

Besides  the  broadening,  one  should  also  account  for  positions  of  the  quantized  energy 
levels  sn  and  for  the  inter-level  spacing  A .  Both  of  them,  sn  and  A ,  depend  on  the  length  of  the 

SWCNT  sections.  The  relevant  frequency  range  corresponds  to  the  THz  domain  (  A/zf;r  =  0.5- 
50  THz,  or,  in  terms  of  wavelength,  AEF  =  1  inm-5  //  m).  This  means  that  the  EF  wavelength 


and  the  nanotube  section  size  match  each  other  just  for  shortest  waves  of  the  THz  domain  only. 
The  range  A^EF  is  limited  by  circuit  parameters  and  by  coupling  to  the  external  environment. 
Besides,  there  are  additional  restrictions  having  an  intrinsic  origin,  as  discussed  in  fonner 
section.  The  upper  frequency  limit  (i.e.  shortest  wavelengths)  is  detennined  by  maximal  level 
energy  difference  A .  The  lower  frequency  limit  depends  upon  coupling  between  the  nanotube 
shoulders  and  the  broadening  due  to  electron-phonon  interaction.  The  inter-electrode  coupling 
depends  on  several  factors  including  transparency  and  geometry  of  the  Schottky  barriers,  and 
also  purity  of  the  normal  electrodes.  Since  the  elastic  electron  scattering  on  the  atomic  defects 
and  on  the  adsorbed  molecules  (which  stick  to  the  nano  tube  sections  [2 -6])  contributes  into  the 
localized  level  width,  it  also  affects  the  lower  frequency  limit  (i.e.,  the  frequency  resolution  of 
the  nanosensor). 

Our  work  suggests  that  the  carbon  nanotube  sensors  have  numerous  advantages  as 
compared  to  other  types  of  devices.  One  important  benefit  is  that  the  intrinsic  properties  of  the 
nanotube  sections  are  controlled  by  gate  and  source-drain  bias  voltages  [5, 6, 8, 9].  In  this  way, 
one  alters  the  electron  bandstructure  and/or  the  charge  carrier  concentration.  The  electron 
bandstructure  is  affected  by  a  transversal  electric  field  [A-D]  in  a  split  gate  configuration. 
Typical  magnitudes  of  the  electric  field  may  be  quite  low  [  E  >  10  V/  p,  m  or  1  mV/A  for  a 
particular  class  of  ’helical'  (n,l)  nanotubes],  or  0.01-1  V  for  other  types  of  semiconducting 
nanotubes  [5, 6, 8, 9].  The  charge  carrier  concentration  and  frequencies  of  quantum  mechanical 
resonances  are  changed  under  the  influence  of  a  bias  gate  voltage  [5, 6, 8, 9].  Using  the  d.c.  bias 
voltages  USD  and  VG ,  one  insures  perfect  tuning  of  the  carbon  nanotube  sensors  as  well  as  a 
prompt  optimization  of  their  working  characteristics.  Another  advantage  of  NJs  is  presence  of 
sharp  intrinsic  resonances  in  their  transport  characteristics.  For  such  reasons  the  carbon  nanotube 
sensors  may  be  easily  tunable  and  work  with  a  good  precision. 

2.  Probing  the  intrinsic  state  of  a  one -dimensional  quantum  well  with  a  photon-assisted 
tunneling. 

The  THz  diapason  of  the  carbon  nanotube  spectral  analyzers  can  be  extended  by  fabricating  of 
combs  the  metallic  carbon  nanotubes.  The  metallic  nanotube  combs  sould  be  combined  with  the 
semiconducting  nanotube  combs  as  shown  schematically  in  Fig.  1.  The  metallic  nanotube  combs 
are  furnished  of  single  wall  carbon  nanotubes  laid  on  a  dielectric  substrate.  Each  of  the 
individual  junctions  in  the  comb  consists  of  the  single  wall  metallic  nanotube  7  with  emitter 
e  and  collector  C  electrodes  attached  to  its  ends  (see  sketch  of  a  single  junction  in  Fig.  7).  The 
junction  infact  furnishes  a  quantum  well  with  quantized  energy  levels  fonned  inside.  Depending 
on  the  tube  properties,  that  setup  (see  Fig.  7)  corresponds  to  various  condensed  matter  systems. 
The  single  wall  carbon  [n,m]  tube  sketched  in  Fig.  5  is  a  rolled  up  atomic  honeycomb 
monolayer  formed  by  two  sublattices  A  and  B.  The  integer  indexes  n  and  m  ( n  >  rn  >  0 ) 
of  the  rollup  vector  R  =  nR,  +  mR.2  actually  determine  the  electronic  bandstructure  of  the 
tube.  In  particular,  if  n  —  m  —  3u  ( v  being  an  integer)  the  tube  is  metallic  while  it  is 
semiconducting  or  insulating  otherwise  [11].  The  charge  carriers  in  metallic  tubes  conform  to  the 
linear  dispersion  law  £k  =  ±t’F  |  k  |  [where  "  ±  "  corresponds  to  electrons  (holes),  and  vF  is 

the  Fermi  velocity].  A  lot  of  discussions  address  the  intrinsic  state  of  metallic  tubes  where  the 
Tomonaga  Luttinger  liquid  state  [14]  (TLL)  may  presumably  occur  [12,13,15,16,17].  In  contrast 
to  semiconducting  tubes,  where  a  general  consensus  is  achieved,  unconventional  features  of  the 
metallic  tubes  are  not  well  understood  yet.  A  lot  of  attention  [12,13,15,16,17]  is  paid  to  the 
strong  correlation  effects,  and  to  the  one-dimensional  transport  of  the  electric  charge  carriers. 
Along  with  the  TLL  state  [14]  in  metallic  tubes  [12,13,15,16,17]  under  the  current  elaboration 
there  are  models  operating  with  non-interacting  electrons  while  other  models  exploit  coupling  of 
the  tube  to  the  external  environment[18].  Although  there  are  indications  of  the  TLL  state  in  the 
shot  noise[15]  and  in  angle  integrated  photoemission  measurements [16],  present  experimental 
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evidences  are  still  indirect[  18,19].  Therefore  more  efforts  to  clearly  identify  the  intrinsic  state  of 
the  one-dimensional  quantum  wells  formed  of  metallic  carbon  tubes  are  required. 


Fig.  7.  (a)  A  quantum 
well  (QW)  composed 
of  the  ID  section  X 
with  emitter  £  and 
collector  £  electrodes 
attached.  The  potential 
barriers  are  shown  in 
black  at  the  <£/X  and 
X/<£  interfaces,  (b) 
Energy  diagram  of  the 
PAT  process  in  the 
QW.  (c)  The  split  gate 
configuration  of  the 
QW.  The  right  side 
inset  shows  how  the 
electric  field  is  applied 
to  the  X  section. 


Typical  quantum  well  setup  [5,6,B-D]  is  sketched  in  Fig.  7  where  the  ID  section  is 
denoted  as  7  .  The  bias  voltage  Vec  drops  between  the  emitter  ( e )  and  collector  ( C  )  electrodes, 
while  the  gate  voltage  VG  is  applied  to  the  n  ++Si  substrate  as  shown  in  Fig.  7.  The 
electrochemical  potentials  in  e ,  7  ,  and  C  are  denoted  as  /i  T  r .  The  e  and  C  electrodes  are 

separated  from  the  metallic  tube  section  7  by  the  interface  barriers  I  shown  in  black  in  Fig.  7 
(a,b).  The  potential  barriers  emerge  from  differences  between  the  Fermi  velocities  in  the  adjacent 
electrodes.  In  this  part  of  work  we  suggest  a  method  which  identifies  the  quantized  levels  of 
charge  and  spin  excitations  in  the  Tomonaga  Luttinger  liquid  state  inside  the  ID  quantum  well 
shown  in  Fig.  7(a-c).  Our  method  exploits  the  fact  that  the  electromagnetic  field  (EF)  interacts 
with  the  charge  excitations  only,  while  the  neutral  particles  remain  unaffected.  When  a  tunneling 
electron  with  energy  e  absorbs  n  photons  of  the  external  electromagnetic  field,  the  intrinsic 
structure  of  the  TLL  state  in  7  is  pronounced  in  the  multiphoton  tunneling  process  probability. 
The  photon-assisted  tunneling  (PAT)  influences  probability  of  the  single  electron  tunneling 
(SET)  which  helps  to  elucidate  the  intrinsic  state  of  the  tube.  In  this  paper  we  address  quantum 
wells  with  long  and  short  7  sections.  When  the  tube  is  long,  the  interlevel  spacing 
A  =  hvF  /  L  ( vF  is  the  Fermi  velocity,  L  is  the  tube  length)  is  small.  Therefore  quantization 
of  the  electron  motion  inside  7  is  negligible.  In  that  case  the  local  single  electron  density  of 
states  A f(e)  inside  7  has  a  dip  at  zero  energy  e  =  0  .  We  will  see  that  such  a  dip  is  clearly 
visible  in  the  photon-assisted  and  in  the  single  electron  tunneling  characteristics  which  helps  to 
identify  the  TLL  state.  In  the  opposite  limit  when  the  tube  is  short,  the  ballistic  motion  of  the 
charged  and  neutral  excitations  inside  7  is  quantized.  During  the  tunneling,  an  electron  splits 
into  four  p±:c t±  -bosons  (two  density  and  two  spin).  The  bosons  populate  quantized  levels  with 
different  energies  ep+  ^  ep  {rjG)  . 


10 


The  charge  boson  energy  levels  are  detected  with  the  photon-assisted  tunneling  (PAT,  see 
Ref.  [B]).  We  will  see  that  the  tunneling  mechanism  is  sensitive  to  the  emitter-collector  Vec 

and  the  gate  VG  voltages.  Therefore  the  TLL  properties  are  pronounced  in  the  differential 
conductivity  curves  of  the  one-dimensional  quantum  wells.  In  measurements  made  on  the  same 
setup,  the  spin  boson  levels  are  fingered  from  Zeeman  splitting  oc  / /B II  /  e  when  applying  a 
finite  d.c.  magnetic  field  H  ^  0.  The  quantization  of  both  the  charge  and  spin  excitations  is 
proclaimed  in  the  differential  conductivity  curves  of  the  quantum  well. 

The  TLL  tunneling  density  of  states  of  a  long  QW 

If  the  metallic  tube  section  7  is  long,  L  >  vftt  (where  rT  =  h  /  Tn  is  the  net  tunneling 
time),  the  quantization  inside  7  is  negligible.  The  level  separation  for  typical  carbon  tube 
junctions  [5,6,8-10]  becomes  indistinguishable  when  L  >  3  //  m.  Strong  electron  correlations 

drive  the  electron  system  into  the  Tomonaga  Luttinger  liquid  (TLL)  state[  12-19]. 

A  simple  evaluation  from  the  bandstructure  calculations  gives  j3G  —  0.005  —  0.03 

depending  on  directions  of  the  rollup  vector  R  .  The  M  ( E )  shape  is  also  controlled  with  LSG 

utilizing  the  split  gate  configuration  [A-C]  as  shown  in  Fig.  7(c).  The  electric  field  in  that  setup 
is  perpendicular  to  the  tube  axis  as  shown  in  the  right  inset  to  Fig.  7(c).  The  split-gate  setup 
allows  driving  the  tube  electron  state  from  semiconducting  to  the  metallic  one.  The  transversal 
electric  field  induces  a  finite  dipole  momentum  directed  in  perpendicular  to  the  tube,  which 
renormalizes  g  as  well. 


Fig.  8.  (a)  (a)  Splitting  of  the  zero-bias  TLL  dip  in  the  differential  tunneling  conductivity 
g(Vcc)  [in  units  of  c‘  vF  N(0)]  in  a  long  CNT  junction  due  to  the  photon-assisted  tunneling. 
Spacing  between  the  zero  dip  and  adjacent  satellite  dips  is  MU  e.  (b)The  Coulomb  staircase 
in  the  PASET  current-voltage  characteristics  I(Vdc)  versus  reduced  voltage  Vdc  (see  text)  of 
the  CNT  junction  in  the  TLL  state  with  g=0.4  under  influence  of  the  a.c.  bias  field  with  the 
amplitude  eV(1  -3.4  (in  units  of  8  =  e2/C)  and  for  a  symmetric  junction  (q  =  0.5).  Curve  A 
corresponds  to  fi  =  4.7,  curve  B  to  Q  =  3.7  and  curve  C  to  Q  =  6.7. 

The  corresponding  alteration  of  the  Luttinger  parameter  g  is  evaluated  using,  e.g.,  results  of 
Refs.  [27-29].  It  gives  g  — >  g  ■  J^1  (VSGd,  /  hvF  J  .  For  a  narrow-gap  semiconducting  7  (see 
Ref.  [30])  and  for  typical  parameters  of  the  tube  quantum  well,  the  split  gate  induced  change  is 
g  — >  g  +  /3sgVsg  where  f3SG  =  0.01  —  0.05  for  different  rollup  vectors.  If  the  transversal 
electric  field  I  SG  /  d  inside  7  is  sufficiently  strong,  one  induces  a  semiconducting-metal 
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transition  [27-29].  The  electronic  properties  of  the  tube  then  switch  from  a  one-dimensional 
narrow  gap  semiconductor  to  the  TLL.  The  time-averaged  conductance  cr(Lec)  of  the  long 
CNT  junction  exposed  to  an  external  electromagnetic  field  is  computed  using  Maxwell  equations 
combined  with  the  Green  function  approach.  We  compute  cr(Vec)  and  the  electric  current 

I(Vdc)  [  Vdc  —  (Vec  —  Vt)p  is  the  reduced  voltage,  Vf  —  (E  —  eF)  /  erj  is  the  SET  threshold 
voltage]  for  the  two  cases  of  interest.  One  limit  corresponds  to  6  «  T  when  the  single¬ 
electron  tunneling  is  not  essential  [  8  =  e2  /  C ,  C  =  Ce  +  C  is  the  net  capacitance  of  the 
double-barrier  junction,  Ce(c)  's  the  emitter  (collector)  capacitance].  The  current-voltage 

characteristics  /(Vdc )  [where  Vdc  =  (Tec  —  V] )//  ]  in  condition  of  the  photon-assisted  single¬ 
electron  tunneling  (PASET)  across  the  quantum  well  in  the  Tomonaga  Luttinger  liquid  state  are 
shown  in  Fig.  8(b).  According  to  Ref.  [23],  the  equilibrium  shape  of  the  /(Fdc)  curves  (quoted 
as  Coulomb  staircase)  is  extremely  sensitive  to  the  double-barrier  junction’s  parameters  such  as 
barrier  transparencies,  capacitance,  symmetry,  purity  of  the  carbon  tube  section,  and  the  energy 
level  spacing.  The  photon-assisted  tunneling  induced  by  the  external  electromagnetic  field 
introduces  additional  features  in  those  curves.  We  have  computed  PASET  curves  for  a  QW  with 
a  long  7  -section  where  the  single  electron  tunneling  takes  place.  The  external  x  -polarized 

electromagnetic  field  induces  an  a.c.  bias  voltage  across  the  junction  as  V ^  cos  Clt .  The  most 
remarkable  elements  of  the  7(F^c)  curves  A-C  in  Fig.  8(b)  are  local  dips  which  originate  from 
an  interference  between  the  zero-energy  TLL  anomaly  pronounced  in  equilibrium  at  e  =  0 
[see  Fig.  8(a)]  and  the  photon-assisted  single  electron  tunneling  processes.  The  Coulomb 
staircase  curve  A  in  Fig.  8(b)  corresponds  to  Cl  =  4.7  ,  curve  B  to  Cl  =  3.7  and  curve  C  to 
Cl  =  6.7  computed  for  g  =  0.4  . 

Identifying  of  the  charge  and  the  spin  boson  energy  levels 

In  a  opposite  limit  when  the  7  -section  is  short,  the  quantized  energy  levels  are  well  resolved, 
since  the  condition  Fe,c  «  A  is  observed  (  Fe,c  are  the  e  <F>  T  and  T  <=>  C  electron 
tunneling  rates,  A  =  hvF  /  L  is  the  interlevel  spacing  inside  7  ).  In  this  Section  we  neglect 
by  the  single-electron  tunneling  (SET)  contribution  [23]  (Coulomb  blockade  phenomena).  That 

is  justified  when  the  temperature  T  is  not  too  low,  T  »  Fe,c  .  In  that  limit  we  use  equation  for 
electric  current  again  but  with  a  different  spectral  density  function  Im/C(e)  which  now  acquires 
a  comb-like  shape.  Due  to  the  spin-charge  separation  in  the  Tomonaga  Luttinger  liquid  (TLL) 
there  are  two  sets  of  quantized  energy  levels  in  a  low-transparent  quantum  well  with  a  short  7  - 
-section.  For  the  QW  transparency  T  =  0.3  (where  T  —  4:nF nLn  /  (hvF) ,  Ln  —  Le  +  Lc, 

Le  and  Lr  are  the  e  and  C  thicknesses  respectively,  vF  =  8.1  •  105  m/s)  one  gets  Fn  pa  0.3 
meV.  For  the  tube  length  L  =  3  //  m  one  obtains  spacing  between  the  quantized  levels  as 

A  =  1  meV.  The  photon-assisted  processes  cause  an  additional  splitting  ~0.6  meV  which 
corresponds  to  the  a.c.  bias  frequency  Cl  pa  1  THz.  Following  to  Refs.  [12,22],  one  defines  the 

transmission  coefficient  as  T(E )  =|  ihGR  ( L,E )  |2 .  We  assume  that  coupling  of  the  single  wall 

tube  segment  7  to  the  external  e  and  C  electrodes  is  weak.  In  this  approximation  we  compute 
the  local  electron  density  of  states  Af(e)  implementing  boundary  conditions  [12]  for  the 
electron  wavefunction  inside  a  short  carbon  tube  section  7  .  Then  the  quantized  energy  levels 
are  firmly  separated  from  each  other  and  resolved.  The  charge  p+  bosons  populate  the  energy 

levels  £%+  =  hvFn  /  Lg  =  nA  /  g  (where  n  is  integer  number),  while  three  other  neutral  p 
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and  cr±  -boson  energy  levels  have  conventional  values  c7/'~'rJ±)  =  hvFn  /  L  =  n  ■  A  .  In  the 
steady  state,  when  the  EF  is  off  (i.e.,  oec  =0),  during  the  tunneling  say,  from  e  to  7  ,  an 
electron  splits  into  four  bosons  as  e  — >■  p+  +  p_  +  u+  +  u_ ,  which  assumes  the  energy 

conservation  as  E  +  r\eV^c  —  e£+  +  e^~  +  e°+  +  —  n(3  +  1  /  g)  A  (  n  being  the 

integer  number).  That  corresponds  to  a  resonance  tunneling  through  the  quantized  TLL  states 
tuned  by  Vec  .  However,  if  the  electromagnetic  field  (EF)  is  on  (  ae  c  0  ),  the  resonance 

tunneling  condition  changes.  That  happens  because  the  a.c.  field  acts  on  the  charge  p+  -bosons 
only,  which  absorb  the  EF  photons  during  the  photon-assisted  tunneling  processes. 


Fig.  9.  Splitting  of  quantized  levels  due  to  the  photon  assisted  tunneling  and  Zeeman 
effect  as  pronounced  in  the  single  electron  density  of  states  N(s)  of  a  short  CNT 
junction,  (a)  A  free  electron  quantized  level  (solid  curve)  for  which  the  charge  p+ 
and  spin  a+  bosons  coincide  splits  by  an  a.c.  electromagnetic  field  (with  spacing 
ccfiQ/  e)  and  by  a  d.c.  magnetic  field  (with  spacing  oc  jUbH  /  e )  simultaneously,  (b) 

The  quantized  levels  of  spin  a+  and  charge  p+  bosons  have  different  energies 
A’ 1  ^  san  +  in  the  TEE  state.  The  a.c.  bias  splits  the  charge  boson  levels  (dashed  curve 
on  the  right)  while  the  d.c.  magnetic  field  splits  the  spin  boson  levels  (dotted  curve 
on  the  left)  only.  The  charge  boson  localized  energy  level  £fn’+  splits  in  the  two 
satellite  peaks  with  spacing  2tiQ  /  e .  Although  the  a.c.  field  has  no  influence  to  the 
neutral  spin  bosons,  the  spin  level  £°+  splits  [12]  in  two  sublevels  spaced  with 
Az  oc  jubH  /  e  (both  in  units  of  A)  due  to  the  Zeeman  effect  when  a  d.c.  magnetic 
field  H  ^  0  is  applied. 

The  photons  do  not  excite  the  neutral  p  and  cr±  bosons  since  they  do  not  interact  with  the 
external  electromagnetic  field.  The  a.c.  field-modified  resonance  condition  depends  on  both  Vec 
and  Q  simultaneously 

E  +  r]eVec  =  +  mMfj  +  e^~  +  e°+  +  =n(3  +  l/g,)A  +  mhQ  where  n  and  m 

are  integer  numbers.  The  external  electromagnetic  field  splits  the  conductivity  peaks  selectively. 
Because  Vec  and  Q  are  bound  by  the  condition  (aceson),  this  imposes  a  constrain  on  the  net 
photon-assisted  tunneling  (PAT)  resonant  current  through  the  quantum  well.  We  extract  g  and 
A  from  the  d.c.  PAT  current -voltage  characteristics.  The  splitting  of  the  charged  p+  boson 
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peaks  by  the  a.c.  field  helps  to  identifying  of  the  Tomonaga  Luttinger  liquid  state.  The  method  is 
illustrated  further  in  Fig.  9  where  we  show  a  single  peak  in  Af(e)  corresponding  to  a  quantized 
free  electron  energy  level  [see  Fig.  9(a)].  For  non-interacting  electrons  ( g  —  1)  the  same  single 
level  splits  either  by  an  a.c.  field  due  to  the  photon-assisted  tunneling  phenomena  with  spacing 
oc  mhQ  (  m  being  integer)  or  by  a  d.c.  magnetic  field  with  the  Zeeman  spacing  oc  /_ibH  The 

situation  is  remarkably  different  in  the  Luttinger  liquid  state  when  g  ^  1  and  the  charge  f>_ 

and  spin  a,  levels  have  distinct  energies  e^+  ^  £°+ .  Then  one  easily  identifies  the  charge 

and  spin  levels  merely  by  applying  the  a.c.  field  and  d.c.  magnetic  field  to  the  same  quantum 
well.  If  a  level  splits  with  spacing  cx  mftfL  by  the  a.c.  field  only  (showing  no  response  to  the 
d.c.  field)  then  it  certainly  is  a  p+  charge  boson  level  ( g  ^  1).  If  it  splits  by  the  d.c.  magnetic 

field  [12]  with  the  Zeeman  spacing  oc  pBH  showing  no  response  to  the  a.c.  field,  then  it  must 
be  associated  with  the  spin  bosons  a+  .  However  if  the  both  a.c.  and  d.c.  magnetic  fields  split 
the  same  level,  then  the  level  belongs  to  the  non-interactive  electrons  ( g  =  1 )  as  had  been  said 
above.  In  this  way  one  perceives  the  charge  and  spin  bosons  in  experiments  when  applying  a.c. 
electromagnetic  field  in  combination  with  the  d.c.  magnetic  field  to  a  carbon  tube  junction.  An 
important  requirement  to  the  experimental  metallic  carbon  tube  quantum  well  samples  is  that 
they  must  be  clean. 


in  units  of  A/e].  (b)The  Coulomb  staircase  in  the  PASET  current-voltage  characteristics 
I(Vdc)  versus  reduced  voltage  VdC  (see  text)  of  the  CNT  junction  in  the  TLL  state  with 
g=0.4  under  influence  of  the  a.c.  bias  field  with  the  amplitude  eV(1)=3.4  (in  units  of  8  = 
e  /C)  and  for  a  symmetric  junction  (r|  =  0.5).  Curve  A  corresponds  to  Q  =  4.7,  curve  B 
to  Q  =  3.7  and  curve  C  to  Q  =  6.7. 

An  electron-impurity  scattering  in  real  samples  leads  to  a  formation  of  additional  pairs  of  combs 
with  different  periods.  Then,  an  identification  of  the  TLL  state  becomes  possible  with  a  mere 
generalization  of  the  method  described  above.  Ratio  of  the  noise  power  to  the  mean  current 
(Fano  factor)  is  computed  as  F  —  ^2nTn(l  —  T  )  /  ^2nTn  where  the  summation  is  performed 
over  the  conducting  channels.  Remarkably,  the  multiphoton  absorption  is  pronounced  in  the 
noise  spectra  as  well.  In  this  way  we  suggest  a  method  of  the  noise  spectroscopy  for  studying  of 
the  photon-assisted  tunneling  into  the  TLL  state. 

Phenomena  considered  in  this  part  of  work  originate  from  a  specific  physics  of  the  charge 
and  spin  carriers,  behaving  like  a  blend  of  four  non-interacting  bosons.  The  Tomonaga  Luttinger 
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liquid  state  occurs  inside  the  one-dimensional  quantum  well  formed  by  a  metallic  single  wall 
carbon  tube.  The  TLL  state  is  tested  with  applying  of  an  external  a.c.  electromagnetic  field  and 
of  a  d.c.  magnetic  field  simultaneously.  The  a.c.  field  splits  the  charge  boson  energy  levels  due  to 
the  photon-assisted  tunneling  while  the  d.c.  magnetic  field  splits  the  spin  boson  levels  due  to  the 
Zeeman  effect.  That  allows  a  mere  identification  of  the  quantized  energy  levels  associated  with 
the  charge  and  spin  bosons  forming  the  TLL  state  in  relevant  experiments.  Besides,  one  also 
detennines  the  quantized  level  spacing  A  and  the  TLL  parameter  g .  The  unconventional 
electronic  and  photonic  properties  of  the  metallic  carbon  tube  quantum  well  can  be  utilized  in 
various  nanodevice  applications  including  THz  field  sensors  and  nanoemitters. 

3.  Directional  photoelectric  current  across  the  bilayer  graphene  junction. 

Electromagnetic  properties  of  the  bilayer  graphene  [31-34]  offer  enormous  opportunities  for 
scientific  research  and  various  nanoelectronic  applications.  They  emerge  in  spectroscopy  of 
bound  and  scattering  states,  in  the  photon-assisted  chiral  tunneling  and  in  direct  probing  of 
strong  correlation  effects.  Potential  applications  include  electromagnetic  field  (EF)  spectral 
analyzers,  receivers,  detectors,  and  sensors[B].  The  crystal  lattice  of  the  bilayer  graphene  [31-34] 
consists  of  four  equivalent  sublattices  of  carbon  atoms  while  the  charge  carriers  behave  there  as 
massive  "chiral  fermions" [3 1,34].  The  chiral  fermions  (CF)  in  bilayer  graphene  have  a  finite 
mass  me  h  ,  like  conventional  electrons  ( e  )  and  holes  ( h  )  in  metals  and  semiconductors  [3 1  -34] . 

The  chirality  relates  the  particles  to  certain  sublattice  and  is  responsible  for  various 
unconventional  d.c.  electronic  and  magnetic  properties  of  the  bilayer  graphene  [31].  In  contrast 
to  an  ordinary  tunneling  through  a  conventional  potential  barrier,  during  the  chiral  tunneling 
(CT)  an  incoming  electron  is  converted  into  a  hole  moving  inside  the  graphene  barrier  in  a 
reverse  direction  as  indicated  in  Fig.  11a  (Klein  paradox[38,39]).  This  yields  a  finite 
transparency  T  ^  0  for  incident  electrons  with  energies  E  below  the  barrier  E  <  U0  ( Ufj  is 
the  barrier  height  energy)  occurring  [31]  at  finite  particle  incidence  angles  cf)  ^  0 .  On  the  other 
hand,  the  steady  state  chiral  tunneling  is  blocked  ( T  =  0  )  in  the  longitudinal  direction  (f>  =  0 . 
The  angle-dependent  transparency  makes  the  chiral  tunneling  being  attractive  for  various 
nanoelectronic  applications  [B].  The  potential  barrier  in  graphene  can  either  be  induced  by  the 
gate  voltage  VG  from  a  Si  gate  slab  or  can  be  fonned  by  three  overlapping  graphene  sheets  as 
shown  in  Figs.  9(c,d).  According  to  Ref.  [31],  the  d.c.  gate  voltage  VG  shifts  the  graphene 

barrier  height,  which  controls  the  chiral  tunneling.  That  process  implies  the  wavefunction  phases 
of  electrons  and  holes  being  interconnected  with  each  other  in  the  graphene.  The  phase 
correlations  during  the  chiral  tunneling  can  also  be  directly  tuned  by  applying  of  an  external  a.c. 
field.  Controlling  of  the  electron  wavefunction  phase  by  an  a.c.  field  had  not  been  accomplished 
yet  and  is  the  subject  of  this  paper.  The  electronic  properties  are  described  by  a  spinor 
wavefunction  T  ,  which  components  depend  on  the  angle  6  between  the  electron  momentum 
p  and  the  x  -axis  (see  Fig.  11).  Similar  spinor  description  had  formerly  been  used  for  Dirac 
fermions  [38]  and  for  relativistic  quasiparticles  in  single-layer  graphene  [31,32],  This  part  of 
work  is  devoted  to  electromagnetic  properties  of  a  bilayer  graphene  junction  shown  in  Fig.  11. 
One  may  expect  that  the  differential  tunneling  conductance  cr(0,Lsd)  of  "clean"  samples 
depends  on  the  angle  (f)  between  the  electric  current  j  and  the  x  -direction  (see  Fig.  1 1).  The 
whole  shape  of  a(4>)  versus  the  source-drain  voltage  Vsd  is  very  sensitive  to  properties  of  the 
bilayer  graphene  barrier.  We  begin  with  computing  of  the  steady  state  cr(0,Vsd)  curves  for  a 
graphene  barrier  biased  by  Fsd .  The  steady  state  results  are  then  utilized  for  studying  of  the  a.c. 
properties.  When  an  external  electromagnetic  field  (EF)  is  applied,  it  strongly  affects  the 
directional  diagram  of  cr{(t)^Vs<i) .  In  particular  we  will  see  that  the  external  electromagnetic  field 
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induces  a  finite  conductance  in  the  straightforward  direction  (0  =  0  ),  which  had  been  blocked 
in  the  steady  state. 


Fig.  11.  (a)  Potential  barrier  in  the  bilayer  graphene  controlled  by  the  gate  voltage  Vg 
and  exposed  to  the  external  electromagnetic  field  EF.  The  scattering  states  inside  the 
barrier  originate  from  conversion  of  an  electron  (e)  to  a  hole  (h).  (b)  Two  coupled 
hexagonal  lattices  with  non-equivalent  carbon  atomic  sites  A,  B,  A ,  and  B  in  the 
bottom  and  top  layers  respectively,  (c)  and  (d)  show  two  possible  setups  of  the  bilayer 
graphene  junction.  The  external  a.c.  field  induces  the  directional  photon-assisted 
resonant  tunnelling. 


That  happens  because  the  electromagnetic  field  affects  the  electron-hole  phase  correlations 
inside  the  graphene  barrier  directly.  In  the  steady  state,  when  the  a.c.  field  is  off,  the  electric 
current  is  fully  suppressed  at  Vsd  <  U0  (for  typical  gate  voltage  VG  =  1  V  and  the  Si09 
thickness  d  =  300  mn  one  finds  [32]  UQ  —  2  meV). 


Photon-assisted  chiral  tunneling 

Here  we  examine  influence  of  an  electromagnetic  field  to  chiral  tunneling  and  discuss  the 
intrinsic  noise.  For  studying  of  the  non-stationary  electric  current  across  the  bilayer  graphene 
junction  we  implement  methods  [21,22].  The  graphene  bilayer  is  modelled  as  two  coupled 
gexagonal  lattices  consisting  of  four  non-equivalent  sites  A,  B  and  A,  in  the  bottom  and  top 
layers  respectively  [see  Fig.  11(a)].  For  graphene  junctions  having  finite  dimensions,  the  motion 
of  chiral  fermions  is  quantized.  The  quantization  imposes  additional  constrains  on  the  directional 
tunneling  diagram.  Permitted  values  of  the  angle  (pn  inside  the  graphene  barrier  are  obtained 
from  boundary  conditions  along  the  y  -direction,  so  the  y  -component  of  the  electron 


momentum  p  =  (hk,  hq) 
which  gives 


is  quantized  as  q  —  im  /  IT  (where  W 


=  arctan 


mr  /  (k'eW) 


is  the  barrier  width), 
where 
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e|(l  —  cos 20)  /  2  . 


k'e  =  \]2m  /  h2^ \e-  UQ 


Fig.  12.  A  steady  state  differential  conductivity  a(c)),VsD)  [<t>  is  the  asimuthal  angle, 
Vsd  is  the  source  drain  bias  voltage  in  units  of  the  graphene  barrier  height  eUo] 
through  the  chiral  barrier  created  by  a  bilayer  graphehe  pellet. 


The  external  field  not  only  splits  the  resonances,  but  also  strongly  affects  angular 
dependence  of  the  chiral  tunneling.  That  happens  because  the  a.c.  field  causes  no  influence  to  the 
y  -component  of  the  electron  momentum  q  since  the  graphene  barrier  is  effectively  one¬ 
dimensional.  The  time  dependence  /^(t)  takes  also  place  when  the  a.c.  field  modulates  the 
graphene  barrier  width  as  D  —*■  DQ  +  D1  cos  Qt .  Splitting  of  the  chiral  tunneling  resonances, 
and  the  angular  redistribution  of  the  electric  current  under  the  a.c.  field  influence  is  better 
pronounced  for  a  finite  barrier  height  U{)  0  and  U  —  U0  +  U1  cos  Qt .  From  Eq.  (phi)  one 

can  see  that  4>'  =  0  if  0  =  0.  However,  if  0^0  one  may  observe  spectacular  phenomena. 
In  this  case  an  external  a.c.  field  induces  a  finite  electric  current  for  an  almost  normal  incidence 
0  m  0  ,  which  was  inhibited  when  the  field  was  off.  When0  ss  0,  the  a.c.  field  actually  causes 
additional  photon-assisted  chiral  tunneling  resonances  to  engage.  The  directional  photoelectric 
effect  (DPE)  may  be  realized  in  two  scenarios.  One  scenario  assumes  that  an  electron  beam 
having  a  finite  angular  width  <50  ^  0  enters  the  graphene  barrier  normally.  A  visible  DPE  can 
be  achieved  in  the  setup  shown  in  Fig.  13(b)  where  the  attached  electrodes  1,3  are  made  of  one¬ 
dimensional  conducting  wires.  If  the  wire  is  much  narrower  than  the  width  of  graphene  stripe  ( 
W  «W),  one  may  consider  the  electric  current  as  a  result  of  one-dimensional  propagation  of 
of  electron  along  the  trajectories  under  influence  of  the  bias  voltage.  Such  method  formerly  had 
intensively  been  used  in  numerous  works  devoted  to  point  contact  junctions.  If  the  electric 
current  is  sufficiently  weak,  the  electrons  coming  from  the  wire  into  the  graphene  stripe 
introduce  a  negligible  disturbance  into  the  electron  spectrum  inside  graphene.  The  translational 
invariance  inside  graphene  is  well  preserved.  In  an  experiment  one  may  use  a  STM  tip  for 
imaging  of  the  electron  wavefunction  in  carbon  nanotube  which  show  a  periodic  pattern.  The 
electrode  1  emits  electrons  under  a  small  but  finite  angle  0  (0<<7r,  0^0)  which 

trajectories  are  focused/defocused  by  the  external  electromagnetic  field  as  indicated  in  Fig. 
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13(b). 


Fig.  13.  The  time  averaged  differential  conductance  cr'(t)  [in  units  of  ao= 

(2eVh)7r(2mUo)  "]  of  a  bilayer  graphene  junction  exposed  to  an  external 
electromagnetic  field  which  modulates  the  barrier  height  U(t)  =  Uo  +  Ui  cosQt.  In  Fig. 

1 1(c)  one  may  notice  a  remarkably  strong  DPE  at  Q/Uo  —  1.  This  corresponds  to  curve 
1  in  Fig.  1 1(d)  where  the  peak  spacing  is  determined  by  the  graphene  barrier  length  D. 

The  frequency  dependence  of  the  transparency  is  governed  by  the  directional  photoelectric 
effect.  A  significant  directional  photo-electric  effect  emerges  even  for  a  relatively  long 
wavelength  A  —  1  mm  -  0.01  //  m  (which  corresponds  to  the  THz  domain)  if  the  condition 

\E  —  UQ  =f^|  <<  j^l  is  met.  The  deviation  angle  (f>r  —  arcsin[J|f?|  /  | E  —  U0  T  fl|sin0] 
inside  the  graphene  barrier  considerably  increases  giving  <f>'  »  0.  This  means  that  an  ideal 
transparency  taking  place  in  the  steady  state  at  0^0  is  redistributed  over  the  angle  (f>  after 
the  a.c.  field  is  applied.  The  transparency  peaks  are  actually  shifted  from  finite  angles  0  ^  0  to 
the  nonnal  incidence  angle  0  =  0.  Another  scenario  involves  an  incident  single  electron  which 
enters  the  graphene  barrier  strictly  in  the  normal  direction  (0  =  0)  under  influence  of  a  high 
frequency  THz  wave.  In  that  scenario  an  electron  absorbs  a  THz  photon  having  the  finite  energy 
E()  and  momentum  q  along  the  y-axis.  Then  the  electron  deviation  angle  <50  just  before 
entering  the  barrier  is  small,  8(f)  «  i r  .  For  instance  taking  v  —  30  THz  (which  corresponds  to 
the  photon  energy  En  =  125  x  10~3eV)  one  gets  8(f>  ~  q  /  k  =  2  x  10-3.  The  photoelectric 
effect  is  well  pronounced  for  an  electron  with  energy  Ee  —  2  •  10  3  eV  after  it  gets  inside  the 
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graphene  barrier.  There  if  \E  —  C/q  =F  the  deviation  angle 

<f>'  —  arcsin[^|ff  |  /  |  E  —  U0  =f  O |  sin  0]  increases  considerably,  since  the  photon  energy  is 
pretty  high,  EQ  /  E  —  50,  =  0.1  eV.  Practically  this  means  that  one  must  set  Ml  ~  UQ  to 

get  a  strong  photoelectric  effect.  In  the  above  example  the  last  condition  also  supposes  that  one 
should  use  U0  ~  En  =  125  meV. 
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Fig.  14.  (a)  The  contour  plot  a(VsD,c|)).  (b)  The  3d  plot  of  a(VsD,<j>)  where  the 
additional  angular  dependence  is  induced  by  an  external  a.c.  field  applied. 


Below  we  consider  two  most  important  field  polarizations  along  the  x  and  z  axes  as  shown  in 
Figs.  9(a,b).  The  barrier  transparency  T(E,  0)  is  affected  by  the  a.c.  field  directly  in  either  case. 
In  particular,  the  barrier  shape  is  modulated  by  the  a.c.  field  polarized  along  the  x  -direction  as 
sketched  in  Fig.  11(a),  since  e  — >  e0  +  eL  cos  (fit).  On  other  side,  if  one  applies  an  a.c.  field 

polarized  as  E  =  (0. 0.  E,) ,  it  modulates  the  barrier  height  since  VG  — >  VG  }  +  Vglj  cos(Q/,  ]  [ 


is  the  steady  state  gate  voltage,  VG  )  is  the  a.c.  field  induced  addition,  see  sketch  in  Fig. 
11(b)].  Then  the  a.c.  field  induced  correction  to  the  d.c.  tunneling  current  is 


h  =  2eId£X£ 


St 


2ne  n£+Q-eV  ne-Q-eV 


)■ 


where  the  transmission  amplitude  t\  f ) 


is 


obtained  from  corresponding  non-stationary  boundary  conditions  at  x  =  0  and  x  —  D. 
Physically,  the  directional  photoelectric  effect  (DPE)  comes  from  an  ingenuous  influence  of  the 
external  electromagnetic  field  to  the  electron-hole  phase  correlations  during  the  chiral  tunneling. 
Technically,  modulation  of  the  barrier  height  by  the  a.c.  field  shifts  positions  of  the  sharp  peaks 
in  the  energy-dependent  barrier  transparency  T(e  ±  fi).  Besides,  it  also  modifies  the  overall 

angular  distribution  of  the  electric  current,  so  the  electron-hole  conversions  occur  with  an 
additional  phase  shift.  Numerical  results  for  both  the  cases  are  presented  in  Fig.  14(b). 
Corresponding  plots  for  the  steady  state  differential  conductance  crJV  .)  and  for  the  photon- 


(P 

—t 


assisted  chiral  tunneling  conductance  =  djx  /  dVsd  —  a  —  aQ  both  indicate  the 

angular  redistribution  of  the  photon-assisted  chiral  tunneling  current  across  the  graphene  barrier. 
The  steady  state  conductance  curve  a(j  in  Fig.  14(b)  corresponds  to  U(l  —  2  meV  while  curve 


ax  is  computed  for  VG  =  1  V  and  Q  =  1  THz.  The  DPE  is  well  illustrated  by  the  sharp 
scattering  resonance  taking  place  in  cr,  (V^)  [see  the  crisp  peak  at  the  incidence  angle 
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Fig.  15.  (a)  A  setup  consiststing  of  a  graphene  pellet  G  with  ballistic  transport  of  charge 
carriers.  The  angular  distribution  of  the  electric  current  is  detennined  using  carbon 
nanotube  electrodes  1-6  attached  to  the  pellet.  The  pellet  and  the  nanotubes  are  laid  on 
the  SiC>2  substrate  and  are  controlled  by  the  gate  voltage  VG+Vc,ac  where  the  2nd  tenn  is 
induced  by  the  external  electromagnetic  field  EF.  The  EF  also  causes  an  angular 
redistribution  of  the  time-averaged  electric  current  which  can  be  used  for  detennining 
the  THz  field  amplitude,  frequency,  and  polarization,  (b)  The  time-averaged  differential 
conductivity  between  the  carbon  nanotube  electrodes  2  and  4  versus  the  field  frequency 
Q  [in  units  of  MU  e  ]  for  different  field  amplitudes  (curves  A,  B,  C,  and  D). 

When  the  a.c.  field  is  off,  the  steady  state  tunneling  at  V(i  =  0.5  in  the  straightforward  direction 
is  suppressed  [see  the  corresponding  curve  aQ  (Vsd  j  for  0  =  0].  However,  if  one  applies  the 
a.c.  field  with  frequency  Q  and  E  =  (E r.  0. 0) ,  it  opens  tunneling  channels  in  the 
straightforward  direction  0  =  0  as  is  evident  from  curve  cr1  in  Figs.  12(b).  In  Fig.  14(b)  we 
compare  two  time-averaged  conductance  curves  <j1  (Vsd  )  under  influence  of  the  a.c.  field  with 
two  different  polarizations  along  the  x  (curve  X)  and  2  (curve  Z)  axes  correspondingly.  In 
either  case  the  cr1  ( Vsd  )  curves  show  remarkable  sharp  peaks,  which  position  however  changes 
versus  the  field  polarization.  Although  the  above  results  are  illustrative,  they  focus  solely  on  the 
limit  of  a  weak  electromagnetic  field  U1  «u0. 

Intrinsic  noise  in  the  bilayer  graphene  junction  originates  as  follows.  The  thennal  noise 
comes  from  the  phonons  emitted  in  the  electron-phonon  collisions.  Matrix  element  of  the 

electron-phonon  collisions  according  to  Ref.  [43,  44]  is  M  ,  oc  ^  p  |  M  ( x )  |  p'  ^  cos  ( 4>pp, )  where 
( f>  i  is  the  angle  between  initial  and  final  states.  The  phase  factor  cos  cppp,  plays  quite  a  different 
role  in  the  bilayer  graphene  compared  to  the  single  layer  graphene  [43]  where  it  is  rather 
cos  ( 0  ,  /  2)  instead.  In  the  latter  case,  the  factor  ensures  suppression  of  the  electron-phonon 

and  electron-impurity  collisions  and  the  transport  of  the  change  carriers  remains  ballistic  up  to 
room  temperatures.  In  contrast,  thennal  noise  in  the  bilayer  graphene  devices  is  rather  high  at 
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room  temperatures.  Another  intrinsic  noise  (Poisson  noise)  arises  due  to  the  "Zitterbewegung" 
effect,  which  is  linked  to  a  jittering  motion  of  the  change  carriers  when  electrons  are  randomly 
converted  to  holes  forth  and  back.  That  produces  noise  even  in  zero  temperature  limit.  The  noise 
is  characterized  by  the  Fano  factor  F  =  )Z((T  ( 1  —  T  j  /  ,  where  T  is  the  tunneling 

probability  in  the  n  -th  channel  and  the  summation  is  perfonned  over  all  the  conducting  channels 
(in  our  setup  this  means  just  integration  over  the  electron  propagation  angle  4>  )•  Our  calculations 

show  that  the  Poisson  noise  becomes  extremally  low  at  F^d  >U0. 

In  conclusion  we  computed  the  electric  current  across  the  bilayer  graphene  junction  in 
conditions  when  an  external  electromagnetic  field  is  applied.  We  have  found  that  the  threshold 
absorption  of  the  external  electromagnetic  field  strongly  depends  on  the  a.c.  field  frequency  and 
amplitude.  The  electromagnetic  field  induces  an  ideal  transparency  of  the  graphene  barrier  in  the 
longitudinal  direction,  which  had  been  fully  suppressed  when  the  a.c.  field  was  off.  That 
directional  photoelectric  effect  originates  from  an  angular  redistribution  of  the  whole 
transparency  diagram  since  the  sidebands  at  finite  angles  are  redirected  to  the  normal  incidence. 
An  experimental  observation  of  the  spectacular  directional  optoelectric  phenomena  would 
provide  a  strong  evidence  for  existence  of  the  massive  chiral  fennions  in  the  bilayer  graphene. 
We  emphasize  that  the  threshold  absorption  emerges  purely  from  a  quantum  mechanical  phase 
shift,  and  not  from  an  inelastic  excitation  by  the  a.c.  field.  That  means  no  heating  is  involved 
during  the  absorption.  The  a.c.  current  induced  by  the  electromagnetic  field  across  the  graphene 
junction  has  a  sharp  angular  dependence,  which  potentially  can  be  exploited  in  sensor 
nanodevices  of  the  external  electromagnetic  field.  The  directional  photoelectric  effect  in  the 
double  layer  graphene  junctions  is  a  unique  phenomenon  which  exists  in  that  system  and  had  not 
been  noticed  in  other  systems,  like  junctions  composed  of  single  layer  graphene  or  of  nonnal 
metals.  Most  intriguing  feature  is  the  switch  between  zero  and  finite  conductance  occurring 
without  energy  absorption.  The  phenomena  considered  above  have  a  great  potential  for  various 
nanoelectronic  applications. 
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